نشریة مهندسی عمران فردوسی سال بیست و هفتم» شماره دی ۱۳۹۵ 
مدل‌سازی بالازدگی تونل و شریان‌های حیاتی در زمین‌های مستعد روانگرایی تحت بار لرزه‌ای * 


علیرضا باقریه!٩‏ ایمان لولوتی(؟ امیرحسین باقریه"" 


چکیده در بسیاری از زلزلاهای اعی رآسیب‌ها ی گسترده‌ای به تونل‌های شهری و حارج شهری وار د آمده است. یکی از عوامل بسیار مخرب 
در این زلزلههاء پدیدة روانگرایی عاک می‌باشد. روانگرایی باعث کاهش تنش موثر و درنتیجه کاهش مقاومت برشی می‌شود و این موضوع 
موجب ایجاد تغییر شک‌های بزرگ در نعاک و ساز مدفون د رآن می‌گردد. روالگرایی موجب جوشش حاک» گسترش عرضصی, فرونشست 
سازه و بالازدگی سازه مدفون می‌شود. در بسیاری از زلزلها بالازدگی به‌علت روانگرایی نحاک اطراف سازه به‌عنوان یکی از عوامل مخرب بر 
تونل‌ها و سازه‌های مدفون گزارش شده است. در این پژرهش رفعار توثل دز زمین‌های مستعد روانگرایی مدل‌سازی شلدء اسسشت. از تحلیان 
تقفاضل محدود دویعدی برای مدل‌ساز ی آزمایش سانتریفیوژ استفاده گردیده است. مدل رفتاری ماس )9 لآبرای مدل‌سازی رفتار ماس 
مستعد روانگرایی در نظ رگرفته شدة است. نتایج مدل‌سازی با اندازه‌گیری‌های صورت گرفته د رآزمایش سانتریفیوژ مقایسه گردیده است و 


توانایی مدل برای ارزیابی پدیدة روانگرایی در اطراف تونل‌ها مورد ارزیابی قرا رگرفته است. 


واژه‌های کلیدی شریان‌های حیاتی؛ روانگرایی؛ بالازدگی؛ زلزلة؛ سانتریفیوز. 
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* تاریخ دریافت مقاله ٩۲/۱۲/۱۳‏ و تاریخ پذیرش آن ٩۳/۱۳/۷‏ می‌باشد 
(۱) نويسنده مسئول: استادیار مهندسی عمران. دانشکده عمران و معماری, دانشگاه ملایر. 
(۲) دانش آموخته کارشناسی ارشد مهندسی عمران - سازه» دانشکده عمران و معماری, دانشگاه ملایر. 


(۳) دانشجوی دکتری مهندسی معدن, دانشگاه علوم وتکنولوژی میزوری» رولا» آمریکا. 


۱۳/۸ 


مقد مه 

شریان‌های حیاتی از آن جهت با این نام خوانده شد؛‌اند 
که قطع‌شدن و پا اختلال در آنها زندگی روزمره بشر را؛ 
به‌خصوص در حین زلزلةٌ با مشکلات عمده روبه‌رو 
می‌سازد. شریان حیاتی نامیده می‌شوند. در سال‌های 
اخیر با رشد و توسعة زندگی شهری. سیستم‌های حمل 
و نقل و تمایل به سازهای زیرزمینی. استفاده از تونل‌ها 
و سازهای مدفون رو به افزایش بوده است. 

از جمله زلزلهایی که آسیب‌های گسترده‌ای به 
تونل‌ها, خطوط لوله و سازهای زیرزمینی وارد کرده 
آننت»می تتوان زلدله لومای تا در منال 1111۹۸1 دز 
کالیفرینا و زلزلهٌ تانگشان در سال 2[1۹7۸]در چین را 
می‌توان نام برد. در زلزلٌ ۱۹۹۵ کوب ژاپن خطوط لولة 
انتقال گاز و آب اسیب جدی دیدند و حدود ۲۳۶ مورد 
نشت گاز از خطوط انتقال گاز گزارش شد. که موجب 
آتش سوزی گسترده‌ای شد[34]. آسیب‌های وارد بر 
خحطوط انتقال گاز در زلزلةٌ چی‌چی در سال ۱۹۹۹ در 
تایوان موجب قطع گاز در بسیاری از مراکز صنعتی 
گردید و تا بهره‌برداری مجدد خطوط لوله. ضرری 
اقتضادی ذر خدود ۲۵ میلیسون دلار آمریکا دوریشی 
داشت[5]. 

مصالح دانه‌ای تحت بارهای لرزه‌ای تمایل به 
فشرده‌شد گی دارند. در توده‌های اشباع به‌علت وجود 
سیال حفره‌ای» عدم زه‌کشی خاک و کوتاه بودن زمان 
بارگذاری. از تمایل کاهش حجم توده خاک دانه‌ای 
جلوگیری می‌شود؛ در نتیجه فشار سیال حفره‌ای افزایش 
می‌پابد و اين افزايش فشار سیال حفره‌ای به کاهش 
تنش موثر و خصوصاً در خاک‌های بدون چسبندگی به 
کاهش مقاومت برشی منجر می‌شود [6]. 

سازهای مدفون در خاک‌های روانگرا در حین 
زلزلةٌ ممکن است دچار فرونشست و پا بالازدگی 
شوند[7]. این سازها به‌علت‌های گوناگون از جمله 
حرکت خاک روانگرا به زیر ساز و افزایش فشار آب 
حفره‌ای دچار بالازدگی می‌شوند[8]. زلزلة ۱۹7۶ آلاسکا 
در ایالت متحده آمریکا [9] زلزلة ۱۹7۶ نیگاتا ژاپن که 


حدود ۸۳۷ از منهول‌های سطح شهر دار بالازدگی 
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شدند و این بالازدگی‌ها به‌علت مسدود کردن سطح 
شهر و ایجاد ترافیک. راه‌های دسترسی و امداد را با مانع 
روبه رو ساختند [10]؛ زلزلهة ۳ کوشپرواوکی که 
بیش از ۰ منهول تاارتفاع ۵ متر از سطح زمین 
بیرون آمدند[11] زلرلة کوبه 1211494 زلزلة ۲۰۰۷ 
چاستو[13) زلزلة ۲۰۱۰ شیلی[14] و زلزلة اخیر ۲۰۱۱ 
توکیو [15] از جمله زئزلةهای می‌باشند که تونل‌هاء؛ 
سازهٌ مدفون و منهول‌ها دچار بالازدگی شده اند. 

تو کایدا و همکاران [16] با استفاده از میز لرزه 
نشان دادند هرچه ضخامت لایه روانگرا بپیشتر شود 
قطر لوله و دانسيتهةٌ نسبی خاک در بالازدگی لوله‌های 
فاضلاب با استفاده از میز لرزه پرداختند. آزمایش‌های 
میز لرزه انجام گرفته توسط کوسیکی و همکاران[18] 
نشان داد که بالازدگی و نشست کل بک فضای قرار 
گرفته در خاک روانگرا براساس شکل ۱ می‌تواند 
به‌صورت ترکیبی از اثرات تغییر شکل جانبی(شکل ۱- 
الف) حرکت سیال حفره‌ای(شکل ۱- ب) و تحکیم 
مجدد(شکل اس رخ دهد. 
ایچی و همکاران [19] نشان دادند سرعت بالازدگی در 
ابتدای بارگذاری بیشتر بوده و ذرات ریز دانه در 
سرعت بالازدگی موثر می‌باشند. 

در مسائل ژئوتکنیکی و اندرکنش بین سازه و 
خاک به‌دلیل اهمیت فوق‌العاده تنش‌های گرانشتی و از 
آنجا که تنش گرانشی حقیقی خاک در مدل‌سازی‌ها با 
مقیاس کوچک قابل مدل‌سازی نمی‌باشد. استفاده از 
آزمایش سانتریفیوژ مناسب‌ترین وسیلة برای مدل‌سازی 
فیزیکی می‌باشد. آدالیر و همکاران[20] در دانشگاه 
تر تن کلضیا شا امستفادف از ارسالشی سا نمی زب 
بررسی طرح لرزه‌ای تونل جورج ماسبی و ظرفیست 
روانگرایی خاک اطراف آن و پاسخ‌های ناشی از 
روانگرایی خاک اطراف تونل پرداختند. در سال ۲۰۰۳ 
لینگ و همکاران [21] بالازدگی خطوط لوله در خاک 
مستعد به روانگرایی را بررسی کردند؛ آن‌ها ضمن 
مدل‌سازی آزمایش سانتریفیوژ لوله با قطر بالا به بررسی 


سال بیست و هفتم» شمارة دی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی» امیرحسین باقریه 


روش‌هایی برای کاهش بالازدگی پرداختند. ساساکی و 
تامورا [22] با بهره‌گیری از آزمایش سانتریفیوژ نشان 
دادند که بیشترین جابه‌جایی و بالازدگی سازءه مدفون 
تقریباً در زمان لرزش اتفاق می‌افند. چیان و مدابهوشی 
[7] تأثیرات عمق, قطرء نوع خاک را در بالازدگی ساز؛ 
مدفون با استفاده از آزمایش سانتریفیوژ بررسی کردند. 
بدلیل پیچیدگی مسائل روانگرایی مطالعات عددی 
محدودی در مورد بالازدگی سازه‌ها و شریان‌های 
مدفون انجام شده است؛ از جمله می‌توان به لیو و 
سانگ اشاره کرد[23] که با استفاده از روش المان 
محدود به بررسی لرزه‌ای ایستگاه مترو مدفون در خحاک- 
های اشباع پرداختند. آنها در این پژوهش برای توصیف 
خاک روانگرای اطراف ایستگاه مترو از یک مدل 
پلاستیسیته که توانایی مدل‌سازی رفتار خاک روانگرا در 
حین بار لرزه‌ای را دارا می‌باشد استفاده کردند. آنها برای 
مدل‌سازی مسئله از نرم‌افزار [۱۸-۹۰۵۵۵/۵6-1( ۲۷ 
استفاده کردند. لینگ و همکاران[24] نتایج حاصل از 
مدل‌سازی عددی آزمایش سانتریفیوژ لینگ و همکاران 
در سال ۲۰۰۳ را که انطباق نسبتاً قابل قبولی با نتایج 
آزمایشگاهی داشت. ارائه کردند. آزادی و میر محمد 
حسینی [25] با استفاده از نرم‌افزار ۳1۸6 به مدل‌سازی 
و بررسی رفتار تونل‌های مدفون در خاک‌های دارای 
قابلیت روانگرایی پرداختند؛ در این مدل‌سازی از 
مدل ۳00 برای مدل‌سازی خاک اطراف لوله استفاده 


شده است. آن‌ها نشان دادند افزايش عمق دفن تونل. 
افزايش میرایی سیستم افزایش زوايةٌ اصطکاک داخلی و 
بهحصوص افزایش زاوية اتساع خاک بالازدگی تونل را 
کاهش و افزايش قطر لوله بالازدگی را افزايش می‌دهد. 
سعیدزاده و هاتف[26] با استفاده از نرم‌افزار5ن۱2 
نشان دادند که افزايش دانسیتهٌ نسبی ماسهٌ و وجود رفتار 
اتساعیء بالازدگی را کاهش می‌دهد و نیز لوله‌های با 
قطر بالا دچار بالازدگی بیشتری می‌شوند. 

هدف از این پژوهش ارزیابی عملکرد و تبیین 
نقاط قوت و ضعف مدل ماسهٌ 1180 در بیان رفتار 


خاک روانگرا در اطراف تونل است. با استفاده از روش 
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تفاضل محدود و به‌کارگیری مدل رفتاری ماسة 7۳0 
آزمایش سانتریفیوژ لینگ و همکاران[21] مدل‌سازی 
شدة و نتایج حاصله با نتایج آزمایشگاهی و نتایج 
شبیه‌سازی المان محدود دوبعدی لینگ و همکاران[24] 


مقایسه شده است. 


مدل‌سازی عددی بالازدگی سازهای زیرزمینی در 
آزمایش سانتریفیوژ 
در این پژوهش با استفاده از نرم‌افزار ۳1 نسخه ۷ به 
شبیه‌سازی آزمایش سانتریفیوژ ارائه شدة توسط لینگ و 
همکاران پرداخته شده است. نرم‌افزار فلک برای 
مدل‌سازی مسائل توأم از تثوری بایوت استفاده می‌کند. 
لینگ و همکاران [21] چندین آزمایش سانتریفیوژ برای 
بررسی بالازدگی و رفتار لرزه‌ای لوله‌های قطور مدفون 
انجام دادند. در این سری آزمايش شتاب گرانش ۳۰8 
اقغال. رهام ایا ماس ی امس نود یا 
دانستیه نسبی تقریبا ۸۳۸ که دارای ذرات با قطر متوسط 
۵ میلیمتر و کمترین و بیشترین وزن واحد حجم به 
ترتیب ۱۳/۸۷ و ۱۷/۳۸کیلونیوتن بر متر مکعب می‌باشد 
و از حلال متولوس با لزجت ۳۰ برابر آب استفاده 
کردند[28]. لوله استفاد شده در این مدل‌سازی از جنس 
آلیاژ آلومینیوم است و دارای طول ۲۸ سانتی‌متر و قطر 
داخلی ٩‏ و قطر خارجی ۱۰ سانتی‌متر می‌باشد(شکل ۲). 
بار دینامیکی در مدل, به‌صورت سینوسی با دامنه ۰/۹۵ 
و فرکانس ۳ هرتز به‌صورت هارمونیک طی زمان ٩‏ ثانیه 
اعمال شد؛ است؛ برای جلوگیری از انعکاس امواج به 
داخل محیط از ظروف لایه‌ای استفاده کردند تا شرایط 
مرزهای نامحدود را فراهم کنند[24]. علی‌رغم آن که 
مدل‌های رفتاری در زمینه‌ی خاک‌هابه خصوص 
خاک‌ه ای روانگرا در حال توسعه و بهبودی 
می‌باشندهیچ یک از آن‌ها نمی‌توانند ادعا کنند که قادر 
به بیان تمام جنبه‌های رفتاری مصالح می‌باشند. هدف 
این پژوهش شبیه‌سازی آزمایش انجام گرفته توسط 
لینگ و همکاران[22] ارزیابی کیفیت مدل ماس 1780 
در پیش‌بینی رفتار یک تونل در خاک روانگرا و تحلیل 
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نقاط ضعف و قوت آن است. به علاوه نتایج مدل‌سازی 
با نتایج مدل‌سازی انجام گرفته توسط لینگ و همکاران 
[24]که با استفاده از روش اجزای محدود و مدل 
رفتاری ۱۷۵۶۲1۲ 2عذهن لمع 02و۲۵ انجام گرفته 
است. مقایسه گردیده است و برتری و کاستی مدل ماسهٌ 
در مقایسه با مدل-1۷21 ۳2۹۲۵۲-71601016۷77102 
111 مقایسه گردیده است. 


معرفی و کالیبره کردن پارامترهای مدل ماسة8] 
در نرم‌افزار برای معرفی پدیده روانگرایی. مدل 
رفتاری 7107 براساس معادلات ارائه شدة بر روی مدل 
استاندار پلاستیک موهر کلمب اشباع که فشار آب 


حفره‌ای براساس کرنش‌های حجمی غیرقابل برگشت 


ب حرکت سیال حفره‌ای 


مدل‌سازی بالازدگی تونل و شریانهای... 


محاسبه می‌شود. گنجانده شده است؛ اما این نرم‌افزار 
قادر به فراخوانی مدل‌های رفتاری علاوه بر مدل‌های 
موجود در خود می‌باشد. در این پژوهش برای 
مدل‌سازی خاک ماسة نوادای اطراف لوله از مدل ماسة 
0 اسفاده شدع است: مدل ماتف 1780 یکا سیل 
پلاسفسته کنش مزثر برای استفاده درانالد های پیش فته 
قرع سکن هش باشستنه فز سفن مات 1736 
پاسخ‌های الاستیک همگن فرض می‌شوند و با مدول 
برشی و مدول بالک به‌صورت زیر تعریف می‌شود [29]: 


6۰-1 
1 (۱) 


۳۳ < 0* (۲) 


ج تحکیم مجدد 


4 ۷۱ ۱۳ ۵ 
و + و۲ سرا ۵ 
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د بالازدگی و يا نشست کلی فضای محدود 


شکل ۱ شکل شماتیک از رفتار سازهُ مدفون در روانگرایی[18] 


یی 3۱ چمتکنام ۱ 
(۱۰۷۲ برمااب ۳‏ مسا 


بته )جیت مزا 


شتاب نگاشت افقی-۸۳1 شتاب نگاشت عمودی ۸۷فشار زمین- لآفشارآب حفره‌ای -۳ 


شکل ۲ شکل شماتیک آزمایش سانتریفیوژ انجام شده توسط لینگ و همکاران [24] 
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۴ و*8 به ترتیسب بیانگر مدول برشی و مدول 
حجمی. 166 ضریب مدول برشی الاستیک. 76 توان مدل 
برشی الاستیک (متغیر بین ۰/4 و ۰0/۱ ۳8 فشار 
اتمسفی 0 بیانگر تتش مثر متوسط و » ضریب ارتباط 
بین مدول برشی الاستیک و مدول بالک الاستیک 
می‌باشد. این ضریب وابسته به ضریب پوآسون الاستیک 
باکت مقلاژ انم ریب هواس هت ۱۳۵ 
یه 

در مدل ماس 17190 کرنش‌های پلاستیک به وسيلة 
سطح تسلیم و قانون جریان کنترل می‌شوند. سطح 
تسلیم در این مدل در فضای تنش, با یک خط گذرا از 
مبدا توصیف می‌شود (مانند مدل موهر کلمب برای 
مصالح بدون چسبندگی). در اولین گام بارگذاری سطح 
تسلیم به وسیلذای مقدار تنش‌های اولیه در آن لحظه 
تعریف می‌شود. مطابق شکل(۳) در حین بارگذاری به 
همان نسبت که تنش برشی افزایش می‌یابد. نسبت 
تنش(۳۲/01) نیز افزايش می‌یابد و منجر به حرکت 
نقطةُ ۸ به نقطة 13 می‌شود؛ در این حالت سطح تسلیم 
نیز تغییر مکان داده و به حط جدید گذرا از نقطه ظ و 
بدا مختصات متفل می‌شود: این تخیر سطح سوم 
منجر به سخت‌شدگی و این سخت‌شدگی توأم با ایجاد 
کرنش‌های پلاستیک می‌باشد. این کرنش‌های پلاستیک 


شامل کرنش‌های برشی و کرنش‌های حجمی می‌باشند. 


٩۵۲ ٩۱۳, ۲ 


۱۵۵۸۵۱ ۱۱۵۰۱۷۵ ۹۱۴69, 


شکل ۳ سطح تسلیم در مدل یوبی‌سی‌سند 


(بیتی و بیرن[29]) 


دارند: 
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۱۳۱ 


۳( 0 
ب 
که *) مدول برشی پلاستیک می‌باشد و با فرض 
یک رابط هذلولی بین 0۷ و 0۳.07 بدین صورت 
تعریف می‌شود: 


2 
| یاه " 
11 


تنش»ع[1 سطح تنش در لحظه‌ی گسیختگی با معرف 
نسبت شکست می‌باشد. 


2 


همان‌طور که گفته شد, کرنش‌های پلاستیک 
شامل کرنش‌های برشی و کرنش‌های حجمی می‌باشند. 
ارتباط بین کرنش‌های حجمی پلاستیک و کرنش‌های 
برشی پلاستیک با استفاده از روش‌های انرژی به‌صورت 
زیرتعریف می‌شود: 


1 ۳ 
ی رع1 
) ۷ 4 0 8 


که ,م۵ زواية اصطکاک داخلی حجم ثابست 
میی‌باشد: 

سطح تسلیم و جهت کرنش‌های پلاستیک 
براساس قانون جریان در شکل () نشان داده شده 
است. در نسبت‌های تنش کم فشردگی پلاستیک رخ 
می‌دهد در حالی که هيچ‌گونه فشردگی در سطح تنش 
متناسب ,0 پیش‌بینی نمی‌شود. برای نسبت تنش‌های 
بزرگ‌تر از ,2 انبساط و یا انقباض برشی پیش‌بینی 
تشه 29 


۱۸ ۳۵۱۵۱۸۵۱ عذاک۳۱۵ 


او , 0۳ 


شکل 4 جهت کرنش‌های پلاستیک همراه با سطح تسلیم 


باربرداري زمانن که مقندار تست کش کاهتن 


سال پیست و هفتم شمار؛ دی ۱۳۹۵ 


۱۳۲ 


می‌یابده تعریف می‌شود و اگر نسبت تنش بدون تغییر 
علامت شروع به افزايش یابد بارگذاری مجدد نامیده 
می‌شود. در بارگذاری مجدد هیچ‌گونه کرنش برشی 
پلاستیک و کرنش حجمی پلاستیک به وجود نمی‌آید و 
تا زمانی که نسبت نش به بیشترین مقدار قبل خود برسد 
بارگذاری مجدد ادامه دارد[29]. 

پارامترهای ورودی مدل ماسةٌ)79] باید براساس 
نتایج آزمایشگاهی کالیبره شوند تا از صحت مقادیر 
پارامترهای ورودی اطمینان حاصل به عمل آید. در این 
پژوهش برای کالیبره کردن پارامترهای ورودی از 
شبیه‌سازی آزمایش برش ساده زه‌کشی نشدة سیکلی در 
نرم‌افزار ۳۲۸ استفاده شدة است. در این شبیه‌سازی» 
تنش قائم برابر با ۱۹۰ کیلوپاسکال و تنش برشی 
به‌صورت افزایشی تا مقدار معینی افزايش می‌یابد؛ سپس 
جهت اعمال تنش برشی تغییر می‌کند و این سیکل تا 
وقوع صفر شدن تنش قانم ادامه دارد. نتایج پیش‌بینی 
روانگرایی با استفاده از آزمایش برش ساده انجام شده 
توسط کانو[30] و شبیه‌سازی عددی در این پژوهش در 
شکل (۵) آمده است. شکل (۵) انطباق مناسب نتایج 
ازهایشگاهی با شیهسازی فد را تشان می ده 
بنابراین می‌توان نتیجه گرفت پارامترهای ورودی کالیبره 
می‌باشند. 

پارامترهای ورودی مدل ماسةٌ)79] باید براساس 
نتایج آزمایشگاهی کالیبره شوند تا از صحت مقادیر 
پارامترهای ورودی اطمینان حاصل به عمل آید. در این 
پژوهش برای کالیبره کردن پارامترهای ورودی از 
شبیه‌سازی آزمایش برش سادة زه‌کشی نشدة نسیکلی قر 
نرم‌افزار ۳۲۸ استفاده شدة است. در این شبیه‌سازی» 
تنش قالئم برابر با ۱۰ کیلوپاسکال و تنش برشی 
به‌صورت افزایشی تا مقدار معینی افزايش می‌یابد؛ سپس 
جهت اعمال تتش برشی تغیبر کرده و این سیکل تا 
وقوع صفر شدن تنش قانم ادامه دارد. نتایج پیش‌بینی 
روانگرایی با استفاده از آزمایش برش ساده انجام شده 
توسط کانو[30]و شبیه‌سازی عددی در این پژوهش در 
شکل (۵) آمده است. شکل (۵) انطباق مناسب نتایج 
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آزمایشگاهی با شبیه‌سازی عددی را نشان می‌دهد. 
بنابراین می‌توان نتیجه گرفت پارامترهای ورودی کالیبره 
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شکل ۵ پیش‌بینی و مقدار پاسخ روانگرایی ماس نوادا 

در شکل‌های (1 و ۷) مسیر تنش مور و رفتار 
تشن کن در یبرع ازشایش سر سانه سرا 
دانسیته‌های نسبی مختلف و نسبت تنش‌های مختلف 


نشان داده شده است. 
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شکل 7 مسیر تنش موّثر(الف) در 121-0.4 


(ب) و نمودار تتش کرنش 
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شکل ۷ مسیر تنش موْثر(الف)در 10.6 (ب) نمودار تنش کرنش 


همان‌گونه که انتظار می‌رود این شکل‌هانشان 
می‌دهند که با افزايش تعداد سیکل‌ها در آزمایش برش 
ساده» کرنش برشی برای رسیدن به یک تنش برشی 
حاص, افزایش می‌یابد و اين بدین معنی است که خاک 
به تدریج با تکرار سیکل‌های بارگذاری و باربرداری: 
مقاومت برشی خود را از دست می‌دهد و در یک تنش 
پرشی مشخص. کرنش برشی بیشتری در آن به وجود 
می‌آید. در قسمت (ب) این شکل‌ها که تنش قائم موثر 
در برابر تنش برشی ترسیم گردیده است. ملاحظه 
می‌گردد که با تکرار و افزایش تعداد سیکل‌های آزمایش 
تنش قانم موثر تدریجاً به‌علت افزایش فشار آب 
حفره‌ای کاهش می‌یابد و در انتهای آزمایش به صفر 
نزدیک می‌شود. این پدیدهُ منطبق بر تعریف روانگرایی 
می‌باشد. همان‌گونه که انتظار می‌رود 038 و 1۳ 


به‌صورت توأم بر تعداد سیکل‌های لازم بر روانگرایی 


سال بیست و هفتم شمار؛ دی ۱۳۹۵ 


۱۳۳ 


شکل (۵) قابل مشاهده است. 


اثر تراکم نسبی ناشی از تتش و معرفی پارامترهای 
خاک اطراف لوله 


همواره دانسیتهٌ نسبی خاک از مهم‌ترین پارامترهای 
معرف ماس می‌باشد. به‌علت تفاوت افزایش تنش در 
قسمت‌های فوقانی و تحتانی بعد از چرخش دستگاه 
سانتریفیوژ تراکم در سطح نمونه کم و در قسمت 
تحتانی مدل بسیار زیاد می‌باشد. پارک و بیرن در سال 
۶ با استفاده از چندین آزمایش تراکم نسبیء نشان 
دادند که در تمام ماسهٌهای مورد مطالعه. تراکم ناشی از 


تنش از رابطه زير پیروی می‌کند[31]. 


0 
م۳ 0+ و2 < ,1 
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که 1۲ تراکم نسبی اولیه و 0۷ تنش موثر 
عمودی و 72 فشار اتمسفر و 0 برابر است با: 


(م.2)01.5-1۲ 9 هگ ان 
۹2 ۵ وه 


من اس هه 


که و6 و «ز6‌به ترتسب بیشترین و کمترین 
نسبت تخلخل می‌باشد و ) عدد سختی ماس است. که 
وابسته به نسبت تخلخل است. این پارامتر براساس 
نتایج پارک و بیرن [31] برای خاک ماسة نوادا (که در 
این مدل‌سازی استفاده شده) ۳۰۰ در نظر گرفته شده 
است. رفتار خاک‌های ماس به شدت وابسته به تراکم 
نسبی خاک می‌باشد با توجه به رابطه 1 دانسیته نسبی 
خاک ماسة وابسته به تنش مور می‌باشد. به‌علت 
تغییرات شدید تنش موژثر در آزمایش سانتریفیوژ و 
وجود ساز؛ با قطر بالا هر المان تقریباً دارای تتش موثر 
منحصر به فردی می‌باشد. بنابراین برای بیان رفتار 
صحیح مسئله. هر المان از مش‌بندی مسئله دارای 


پارامترهای منحصر به فردی برای معرفی مدل 
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اهلاس باشل هی ایس شرهشن با استفاهه از خیتان 
پارامترهای منحصر به فرد وابسته به دانسیتة نسبی آن 


انا اسف آزیه: 


معرفی سیال حفره‌ای 

نرم‌افزار تفاضل محدود ۳1۵0 با استفاده از مشخصات 
سیال موجود در خاک‌های اشباع و برخی پارامترهای 
خاک اشباع با استفاده از تشوری بایوت قادر به 
مدل‌سازی هم‌زمان سیال و خاک می‌باشد. برای در نظر 
گرفتن سیال باید پارامترهای هم‌چون نفوذپذیری خاک 
اشباع. تخلخل خاک. چگالی سیال. درجه اشباع و مدول 
بالک سیال را به‌طور دقیق وارد کرد[32]. 

در اين پژوهش. چگالی سیال را برابر با 18/۳3 
۰ و نسبت تخلخل خاک مطابق با ویژگی خاک 
مورد نظر برابر با ۰/:۲6 در نظر گرفته شده است[27]. 
یکی از پارامترهایی که نقفش مهم و موثری در آنالیز 
مسائل همزمان خاک و سیال دارد. مدول حجمی سیال 


می‌باشد. این پارامتر بدین صورت معرفی می‌شود: 
۸۳ 3 


)۸ ۵۷ 
که ۳ تفییرانت: فعساراور - تست کتریتن 

حجمی می‌باشند. به‌دست آوردن مقادیر مدول حجمی 
با استفاده از رابطهٌ فوق تقریبا غیرممکن می‌باشد [33].در 
بسیاری از آنالیزها و مدل‌سازی‌ها برای معرفی مدول 
بالک سیال از بیشترین مقدار آن یعنی ۲0۳۸ استفاده 
مشود ویبا تام یک مق تجریی قلان کیتیی ز 
ایشا اش مرف مقر خی سیال استشاژه 
می‌شود. به‌علت وجود حباب‌های هوا درون خاک اشباع 
غالبا باید مقدار کمتری از مقدار بيشینه در نظر گرفت. 
کمترین مقداری که برای مدول بالک سیال می‌توان در 


نظر گرفت برابر با مقدار زیر می‌باشد: 
)۹( 0و6 2 1 
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که 0 یک عدد کوچکتر از یک می‌باشد. مقدار 
طول المان و ,0 برابر با چگالی سیال و 8 شتاب گرانش 


وشن 


مرندی و راستی در سال ۲۰۱۲ ضمن مدل‌سازی 
تونل مدفون در خاک روانگرا؛ تأثیر مدول حجمی را 
بررسی کردند. نتایج آن پژوهش تأثیر فوق‌العاده ایین 
پارامتر بر نتایج مدل‌سازی را نشان می‌دهد [34], 

با استفاده از قانون بویل و فرض یکسان بود فشار 
در سیال و هوا. مدول حجمی سیال حفره‌ای تابعی از 
فشار مطق سیال و درجه اشباع. به‌صورت زیر 
می‌باشد [35]: 


> 8,2 
۱۰ و( م1-5) 


۴برابر با فشار مطلق سیال و م5٩‏ درجه اشباع 
اولیه قبل از شروع به چرخش سانتیریفیوژ می‌باشد. 

ماظن که از فتادلی ۱۸ مشیضی اه ها زاین 
مدول حجمی سیال در عمق محیط آزمایش تغییر 
می‌کند و هرچه عمق بیشتر شود. مدول حجمی سیال 
نیز افزايش می‌یابد؛ در نتیجه در صورت یکسان بودن 
کرنش حجمی اضافه فشار حفره‌ای بیشتری تولید 
می‌شود. درجه اشباع یک پارامتر مهم و تعیین‌کننده در 
این تابع می‌باشد. اگر مقدار مدول حجمی سیال بیشتر 
از 7۱۰۳۳۶ ۵ در نظر گرفته شود باید از بیشتربودن 
درجه اشباع از 4۹/۹ اطمینان حاصل شود که چنین 
درجه اشباعی در مسائل عملی غیر واقعی می‌باشد. 
تأثیر پارامتر مسدول حجمی سیال بر بالازدگی در 
شکل(۸) آمده است. با توجه به معادله و تغییرات مدول 
حجمی سیال در عمق. در نرم‌افزار 81۸6 با استفاده از 
زبان برنامه‌نویسی 30 به المان‌های دارای تراز ارتفاعی 
متفاوت مقادیر متفاوتی اعمال شدة است: 


سال بیست و هفتم شماردی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی» امیرحسین باقریه 
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شکل ۸ حساسیت میزان بالازدگی به تغییرات مدول حجمی سیال 


از پارامترهای مهم در مدل‌سازی آزمایش‌های 
سانتریفیوژ و مسائل روانگرایسی نفوذپذیری خاک 
ماسه‌ای می‌باشد. در این مدل‌سازی از ضریب تراوایی 
ماسهٌ نوادا به‌دست آمده در پروژه ولاکس[27] به روش 
آزمایش نفوذپذیری هد ثابت استفاده شدهُ است. با 
توجه به مفاهیم ضریب نفوذپذیری شتاب گرانش و 
غلظت سیال بر ضریب نفوذپذیری تأثیرگذار می‌باشند. 
در آزمایش سانتریفیوژ با لا برابر کردن شتاب گرانش؛ 
ضریب تراوایی لا برابر و با ۷ برابر کردن 
ویسکوزیته سیال ضریب تراوایی آن می شود. 
(۱ 1 ِ 

که 1 ضریب هدایت هیدرولیکی در سانتریفیوژ 
وا ضریب همدایت هیدرولیکی در شتاب 18 
می‌باشد[35] 


مدل‌سازی لوله مدفون (خصوصیات و مشخصات) 
در این پژوهش برای مدل‌سازی لوله آلومینیمی از المان 
لاینر استفاده شدهُ است. المان لاینر توانایی مدل‌سازی 
رفتار الاستوپلاستیک را دارا می‌باشد. این المان برای 
مدل‌سازی المان‌های مصالح غیرخطی مانند دیوار تونل 
و يا مصالح شکل‌پذیر مانند فولاد مناسب می‌باشد. 
پارامترهای المان لاینر به‌صورت جدول (۱) تعریف 


شده است. 


سال بیست و هفتم, شماره دی ۱۳۹۵ 


۱۳۵ 


جدول ۱ پارامترهای لوله استفاده شده در تحلیل عددی 


پارامتر مدل مقدار 


وزن واحد حجم 0۵/۳ | ۲۷/۲ 


مدول یانگ (۱۳2) ۳۱۹۰ 
نسبت یواسون»( ۰/۳۸۹ 
ضخامت لولهه(صه) ۱۵ 


برای مدل‌سازی صحیح رفتار متقابل خاک و لوله, 
از یک المان 
سطح مشترک اجازه جابجای نسبی المانهای خاک را 
نسبت به تونل می‌دهد. در اين پژوهش سختی نرمال و 
پرشی برابر با ۱۲۰ مگاپاسکال بر متر و زاويةٌ اصطکاک 
۵ درجه در اطراف سازهٌ مدفون و خاک در نظر 
گرفته شده است[32]. 


مقايسة نتایج مدل‌سازی عددی لوله‌های مدفون 


در این پژوهش به مدل‌سازی آزمایش سانتریفیوژ با ابعاد 
واقعی پرداخته شدهُ است و برای ارزیابی از صحت 
مدل‌سازی ضمن بررسی بالازدگی سازهُ مدفون؛ اضافه 
فشار آب حفره‌ای در اطرف مدل بررسی شده است. در 
شکل )٩(‏ نتایج حاصل از بالازدگی لوله به روش 
تفاضل محدود بر مبنای مدل ماسهٌ 2136 با مدل‌سازی 
المان محدود انجام گرفته به وسیلهٌ لینگ و همکاران با 
بهره‌گیری از مدل ۷۵/111 ۳5)0۲-21691067162 و 
نتایج آزمایشگاهی مقایسه شده است. شکل )٩(‏ نشان 
می‌دهد مدل‌سازی تفاضل محدود انجام شده در ایین 
پژوهش انطباق مناسبی هم با نتایج آزمایشگاهی و هم 
با نتایج مدل‌سازی المان محدود لینگ و همکاران دارد. 
جابه‌جایی نهایی به‌دست آمده از مدل‌سازی با استفاده از 
روش تفاضل محدود با جابه‌جایی نهایی به‌دست آمده 
از آزمایش سانتریفیوژ انطباق بسیار مناسبی دارد. در 
نتایج این پژوهش, به مانند آزمایش سانتریفیوژ به 
محض اعمال نیروی لرزه‌ای» لوله شروع به بالازدگی 


نشریه مهندسی عمران فردوسی 


۱۳۶ 


می‌کند اما در مدل‌سازی انجام گرفته توسط لینگ و 
همکاران پس از حدود یک ثانیه بعد از لرزش, لوله 
شروع به بالازدگی می‌کند. شیب نمودار مربوط به 
مدل‌سازی پژوهش حاضر نیز با نتایج حاصل از آزمایش 


سانتریفیوژ هم‌خوانی مناسبی دارد. 
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شکل ٩‏ مقادیر بالازدگی سازه مدفون با روش‌های عددی و 
ِ خر 
ازمایشگاهی 


در شکل (۱۰) اضافه فشار آب حفره‌ای در 
پیزومتر ۳1 که در بالای لوله واقع است نمایش داده 
شدة است. نتایج مدل‌سازی این پژوهش در ابتدای 
بارگذاری رزه‌ای انطباق مناسبی با نتایج مدل‌سازی 
لینگ و همکاران دارد اما با مقادیر اندازه‌گیری شده 
پژوهش حاضر به مقادیر اندازه گیری از نظر روند و 
مقادیر نزدیک شدهٌ است. علت اصلی این اخحتلاف در 
لحظات ابتدایی را در مقدار مدول بالک سیال برای زمان 


آغازین بارگذاری لرزه‌ای می‌توان دانست. 
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شکل ۱۰ اضافه فشار آب حفره‌ای در بالای لوله (۳1) با 


روش‌های عددی و آزمایشگاهی 


نشریة مهندسی عمران فردوسی 


مدل‌سازی بالازدگی تونل و شریانهای... 


در شکل‌های(۱۱ و ۱۲) اضافه فشاره‌ای حفره‌ای 
در برابر زمان برای پیزومترهای ۳2 و ۳4 که در کناره 
های لوله واقع شده‌اند. ترسیم ک دیده استه: .در پیزوفتر 
2 شکل (۱۱) ملاحظه می‌گردد که در لحظات آغازین 
بارگذاری لرزه‌ای نتایج مدل‌سازی لینگ و همکاران بر 
مبنای مسدل -۳25/0۲ 27-11 71601067710627 
نوسانات شدیدی را نشان می‌دهد در حالیکه در پاسخ 
های عددی پژوهش حاضر چنین نوساناتی مشاهده 
نمی‌شود و علی‌رغم عدم انطباق مختصری در ابتدای 
بارگذاری از ثانیهةٌ چهارم پاسخ‌های عددی بر مقادیر 
انلازه گیرش شندة انطیاق کامل دارند انم دز ال استق 
که در نتایج مدل‌سازی لینگ و همکاران حتی در 
لحظات آخحرین بارگذاری روند افزايش فشار آب 
حفره‌ای بر خلاف مقادیر اندازه گیری شده ادامه دارد و 
مقادیر فشار آب حفره‌ای نهایی با مقادیر اندازه‌گیری 
شدة تفاوت قابل توجهی دارند. 

در پیزومتسر ۳4 شکل (۱۲) نیز نتایج 
مدل‌سازی‌های عددی در ابتدای بارگذاری با فشارهای 
اندازه گیری شده متفاوت است. مشاهده می‌شود که در 
نتایج پژوهش حاضر نوسانات شدید به‌مانند نتایج 
مدل‌سازی لینگ و همکاران به وجود نیامده است. به 
علاوه فشار آب حفره‌ای نهایی در شبیه‌سازی عددی در 
پژوهش حاضر با مقدار اندازه‌گیری شدهٌ کاملا منطبق 
است. بنابراین می‌توان این گونه جمع بندی نمود که در 
مورد نقاط مجاور لوله شبیه‌سازی با مدل ماس 17130 در 
مقایسه با مدل]11- ۱۷211 ۲25)0۲-216010167102 بهتر 
عمل کرده و فشارهای نهایی در زمان وقوع کامل 
وراک انش یدای بر گتازی تایآ سر مشادیر 
اندازه گیری شدهٌ سانتریفیوژ منطبق است. 


سال بیست و هفتم. شمارفدی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی» امیرحسین باقریه 
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شکل ۱۱ اضافه فشار آب حفره‌ای در سمت چپ لوله (2) با 


روش‌های عددی و آزمایشگاهی 
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شکل ۱۲ اضافه فشار آب حفره‌ای در سمت راست وله (۳4) با 


روش‌های عددی و آزمایشگاهی 


در شکل ۱۳ اضافه فشار آب حفره‌ای در زیر لوله 
نشان داده شده است. اضافه فشار آب حفره‌ای در زیر 
لوله حاصل از مدل ]۱۵/11 ۳29۲۵۲1601060102 
دارای نوسان‌های بسیار زیادی می‌باشد. که با روند نتایج 
آزمایشگاهی انطباق ندارد. اما در مدل‌سازی این 
پژوهش, این نوسان مشاهده نشده است و انطباق خوبی 
در روند و اضافه فشار نهایی بانتایج آزمایشگاهی 
وجود دارد. وجود اختلاف در مقادیر اوليةٌ اضافه فشار 
آب حفره‌ای در زمان‌های اولیه را می‌توان به حساسیت 
مسئله به متغیر بودن مدول حجمی آب به‌خصوص در 
زمان‌های اولیه مرتبط دانست. 


سال بیست و هفتم, شماره دی ۱۳۹۵ 
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شکل ۱۳ اضافه فشار آب حفره‌ای در پایین لوله (3) با 


روش‌های عددی و آزمایشگاهی 


در شکل‌های (۱۶و۱۵) اضافه فشار آب حفره‌ای 
در بالای لوله با نتایج آزمایشگاهی و عددی مقایسه شده 
است که نشان‌دهنده انطباق مناسب نتایج مدل‌سازی 
عددی و نتایج آزمایشگاهی از ابتدای شروع بارگذاری 
لرزه‌ای تا انتهای آن است. در این نقاط نیز مدل‌سازی 
پژوهش حاضر از مدل‌سازی لینگ و همکاران در 


مجموع نسبتا موفق‌تر عمل کرده است. 
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شکل ۱۶ اضافه فشار آب حفره‌ای در فاصله ۱۵ سانتی‌متر از بالای 
لوله در آزمایش سانتریفیوژ (8) با روش‌های عددی و 
آزمایشگاهی 


نشرية مهندسی عمران فردوسی 
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شکل ۱۵ اضافه فشار آب حفره‌ای در فاصله ۱۵ سانتی‌متر از پایین 
لوله در آزمایش سانتریفیوژ(۳9) با روش‌های عددی و 
آزمایشگاهی 


شکل (۱7) تغییرات فشار آب حفره‌ای در 
بالاترین پیزومتر (۳5) را نشان می‌دهد. مشاهده می‌گردد 
که اندازه گیری‌های آزمایش سانتریفیوژ کاهش تدریجی 
فشار آب حفره‌ای در اثر اعمال بار لرزه‌ای را ثبت 
کرده‌اند؛ این روند. در مقادیر حاصل از مدل‌سازی 
عددی پژوهش حاضر مشهود است. بر خلاف این روند 
» نتایج حاصل از مدل‌سازی با روش اجزا محدود 
براساس مدل رفتاری ]11- ۱۷۵71 ۳۵9)0۲-1601016۷102 
در پژوهش لینگ و همکاران[24] نشان دهندهة افزايش 
مختصر اضافه فشار آب حفره‌ای در طی فرآیند 
بارگذاری است. علت اصلی اختلاف بین مدل‌سازی به 
روش المان محدود و نتایج آزمایش‌گاهی عدم 
مدل‌سازی مسائل با کرنش‌های بزرگ و در نتیجه عدم 
پیش‌بینی تورم خاک در سطح زمین در بالای که منجر 
کاهش سطح آب زیرزمینی در مدل‌سازی المان محدود 
می‌شود می‌توان دانست. مدل‌سازی با استفاده از نرم‌افزار 
به روش تفاضل محدود فادر به مدل‌سازی 
کرنش‌های بزرگ می‌باشد و می‌تواند تورم خاک در 


نشریة مهندسی عمران فردوسی 


مدل‌سازی بالازدگی تونل و شریانهای... 
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شکل ۱۰ اضافه فشار آب حفره‌ای در عمق دفن ۱ متری (۳5) 
روش‌های عددی و آزمایشگاهی 
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شکل ۱۷ اضافه فشار آب حفره‌ای در عمق دفن 9 متری (۳6) با 


روش‌های عددی و آزمایشگاهی 


8 10 


در مورد پیزومتر ۳6 که پایین ترین پیزومتر واقع 
در زیر لوله می‌باشد. نتایج مدل‌سازی‌های عددی و 
اندازه‌گیری‌های سانتریفیوژ در شکل (۱۷) به‌صورت 
هم‌زمان نشان داده شده است. در این پیزومتر نتایج 
مدل‌سازی پژوهش از لحظة آغازین بارگذاری لرزه‌ای 
بر مقادیر اندازه گیری مطابقت دارد. در مقابل نتایج 
مدل‌سازی لینگ و همکاران دارای نوسانات زیاد در 
فشار آب حفره‌ای این پیزومتر است. بنابراین می‌توان 
گفت مدل ماس 90[ در اين پیزومتر موفق‌تر عمل 
کرده است. 

در شکل (۱۸) اضافه فشار آب حفره‌ای در پایان 


بارگذاری لرزه‌ای در مقابل عمق ترسیم گردیده است؛ 


سال بیست و هفتم. شمارفدی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی» امیرحسین باقریه 


این شکل مطابقت مناسب مقادیر مدل‌سازی پژوهش 
حاضر با مقادیر اضافه فشار آب حفره‌ای ثبت شده در 
آزمایش سانتریفیوژ را نشان می‌دهد. همان‌طور که 
مشاهده می‌شود اضافه فشار آب حفره‌ای با افزایش 


عمق افزایش می‌یابد. 
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شکل ۱۸ توزیع اضافه فشار آب حفره‌ای در برابر عمق در پایان 


اعمال بار لرزه‌ای 


شکل ۱٩‏ بردار جابه‌جایی المان‌های اطراف لوله 


در شکل (1۹) بردارهای جابه‌جایی المان‌ها در 
پایان آنالیز نشان داده شدة است. همان‌طور که گفته شد. 
یکی از دلایل بالازدگی سازه مدفون» حرکت خاک 
روانگرا به زیر سازهُ می‌باشد. همان‌طور که از شکل قابل 
برداشت می‌باشد. مطابق نتایج آزمایشگاهی حرکت 
خاک روانگرا به سمت زیر سازهٌ مدفون در بالازدگی 


آب حفره‌ای ایجاد شده در حین بار لرزه‌ای دانست. 


سال بیست و هفتم شمار؛ دی ۱۳۹۵ 


۱۳۹ 


نتیجه گیری 

در این تحقیق مدل ماسه 118 برای شبیه‌سازی رفتار 
یک خاک مستعد روانگرایی در اطراف تونل به کار 
گرفته شدة است. نتایج یک آزمایش سانتریفیوژ برای 
ارزیابی صحت عملکرد مدل در بیان رفتار خاک مورد 
استفاده قرار گرفت. در عين حال نتایج با مدل‌سازی 
عددی بر مبنای مدل ]۷۵۲-11 ۳۵۵۱۵۲-216010007102 
مقایسه گرداید. 

مدل‌های رفتاری غالباً بر مبنای آزمایش‌های برش 
ساده و سه محوری ارائه و کالیبره می‌گردند؛ از این‌رو 
اگر یک مدل رفتاری برای شبیه‌سازی چنین آزمایشی به 
کار روت اخبال تطانی مدل‌سازی با مقادی اندازه گر 
شده بیشتر از سایر مسیرهای تنش خواهد بوده این در 
حالی است که در یک محیط پیوسته مانند خاک اطراف 
تونل. مسیرهای تنش با آنچه در آزمایش‌های 
آزمایشگاهی متعارف اعمال می‌گردد. متفاوت است. در 
این پژوهش برای کالیبره کردن پارامترهای مدل ماسة 
6 آزمایش برش ساده انجام شده بر روی خاک 
مورد مطالعه به وسیلة نرم‌افزار ۳1,۸6 شبیه‌سازی 
گردید. 

نشان داده شد که وابستگی تراکم نسبی خاک به 
تنش محدود کننده و نیز مدول حجمی آب به فشار آب 
حفره‌ای بر نتایج مدل‌سازی تأثیر قابل توجهی دارند؛ از 
این‌رو اين پديدةٌ در آنالیزها اعمال گردید. همان‌طور که 
انتظار می‌رود هیچ مدل رفتاری برای زئومتری ادعای 
انطباق مطلق نتایج شبیه‌سازی با آزمایش را در تمام 
مسیرهای تنش ندارد. در مجموع انتظار نمی‌رود مدل 
رفتار در بیان رفتار پیچیده‌ای مانند روانگرایی در تمام 
نقاط و تمام پارامترها کاملاً منطبق بر واقعیت عمل کند. 
به علاوه در مورد اندازه‌گیری‌های سانتریفیوژ نیز 
نمی‌توان ادعا کرد که کاملاً صحیح و عاری از خطاهای 
سینماتیک و ابزاری باشد. هدف از مقايسة مدل‌سازی با 


نشریه مهندسی عمران فردوسی 


۱۴۰ 


آزمایش بررسی قابلیت کلی مدل در بیان رفتار حاک 
است. عواملی مانند مدول حجمی آب (شکل ۸ 
به‌علاوه مرزهای مدل سانتریفیوژ که به‌صورت لایه‌ای 
می‌باشد و ممکن است به شکل ایده‌آل عمل نکند و 
موجب عدم تطابق مدل‌سازی با اندازه‌گیری شود. با این 
توصیف نتایج مدل‌سازی انجام شدهْ همخوانی مناسبی 
با اندازه‌گیری‌های صورت گرفته در آزمایش سانتریفیوژ 
دارند. این موضوع به خصوص در مقایسه نتایج 
مدل‌ساز ی ها با نتایج مدل ۱۷۵ ۳۵۵)۵۲-1601060/1627 
1- انجام گرفته توسط لینگ و همگاران جلب توجه 
می‌کند. بنابراین می‌توان این‌گونه نتیجه گرفت که 
عملکرد مدل ماسه 90] بهتر از مدل مذکور است. 


بر 


مدل‌سازی بالا زدگی تونل و شریانهای... 


آنالیزهای انجام شد؛ د این پژوهش از انیه چهارم 
انطباق مناسبی با نتایج اندازه گیری شد دارند. که علت 
عدم تطابق اولیه عدم معرفی و اطلاعات کافی در مورد 
مدول بالک سیال (بهحصوص اطراف لوله) در لحظات 
ابتدایی می‌باشد. از این‌رو نتایج تحقیق را می‌توان گامی 
جلو برای رسیدن به مدل‌سازی‌های دقیق‌تر در 
خاک‌های مستعد روانگرایی ارزیابی نمود. 

با توجه به نتایج حاصل از اين پژوهش می‌توان 
نتیجه گرفت که در مدل ماسة 11۳80 قابلیت‌های مناسبی 
در مدل‌سازی رفتار خاک مستعد به روانگرایی وجود 


دارد. 


مطا ۶ وامعموه لهمتصطعاممم 24 ماما ۲.۷۷۰ رهقهعظ صه .۲۱.۲ اه0۳ا9 .1 ,000۷0 وم رقتتام 0 .1 


1 5عمصد0ظ مهم چم عمصمتعتصه اهممتاهصعاما 24 که فمصتلععم۶۵ظ "رمعلهنمطانجظ مامتزظ محصمن 1989 


۰ ,۱۷۲۵ رهاامک-ناموو ۵ اهتنا روم‌تصهصونا ارم مه ممتزممصتمم ملمتوطاتهط احمنصدهع‌اوعت0 


1601-1612,)1991( 


ص ملمیامطنه۴ ممطوع‌صه 1 4صه ممهطمنهامصتسن ممتامه]مناونا اامو وی مقامناعاه1 2040 ور ۳06008200۲8 .2 


.(1984) ,11-29 00 ,4 .۱0 ,24 ,۷۵۱ مقطم02ص۲ 8صع علامو ۲ ,5۵۷1۵ ۸ بعصئط 


(10]602010021:)1995 ۲۵۳۲ .معلقناوطانهه وجام 1995 .17 صمیاصول فط1 بت0موو۲ تتقصصصتاه م۲ .3 


60 بتااه از ممم۳۱۷۵۵ ,1995 نامه 17 موم وتفمتتام۳۳۳ وباظ. هصو ۲‏ ین صمطا۹62۳ .4 


٩۱0.3, 00. 146-157, )1995(.‏ و17 ۷۵ : وعتنامناتاگ مهم صمم ۰ ملمنامطاتوم 


ملمنمطانهه 1-1 مط صا ومم‌مصصع موی مصتلومرنم قمع آهسه۳(۲ ۷۲۵۰ ممط ,0۷۷ ۷۷۱۲ ,ظ منطو ,۷۷۷۷ رصع .5 


2000(۰) ,2008008 مصصععو وم گصمی قصمتامصه‌اصا طا6 مط ۵۶ ۳۶۵۵ حص .(بامتصول ۵۶ بدا 16) 


۳۲۵56 1 مطمتانظ ۹۵6080 ردطمهمتوود لها آمعتاتن ۸ :صمتامجکمناوانا (ا90. وکا مطعظ ,1۷ ,101107165 .6 


)2015(: 


مصمومص ۵۶ اواج تمامصصهنل مه طاومل میاه گم )۳۴۲۵۵ تطفیاطماه۷۵0 ٩۳.0.‏ مصقتط ).و .7 


-181 ط ,41 ۷۵۱ ,عصلتهءصنهم ملمنمطاتجط مه ومتصهمرنا لامک . "فلامع. 0علمنام11_ صا. وه‌متتام‌ناتاه 


متصوله٩‏ (1۳1) وبا روطاعمهر1. مممطوگگن اصمصصووعوفه ممتاهم1)رفصه ول راتممرفک ص۳۳1 (2008) متهنا۲ . 


190,)2012( 


[مه 


07 ۱2/0 . 1 خفظ مصمتاجمممن مستامنتعع لصا ماع مه 0متوم۵تر ‏ وتامع[۳۳۵ 0111تا؟۴ 


2008(۰) ۲۱1۷ رحوتع ۲۱۵ 


ممنمطانه۲مصتاهلنا "رومصتاهموتم گم همان مه تماهدهه ‏ متصوزمو" . ,متام 0 28 .ل.۲۱۵1[:۷۷ ,9 


نشریة مهندسی عمران فردوسی 


سال بیست و هفتم شماردی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی, امیرحسین باقریه ۱۳۱ 


(761-773,)1991 .00 ,۸9۳۲ و(:60) 0تفعوهن رعصلهعص1مط۲ 

مج ماما ماهمامه۳ ۵۴ مممصهل ۰۵۵۹۵۲۷۵۵ ۷۲۰ رتطفتومک 20 بن) بل رقتقکا ریگ بلق و. ما10 .10 
۵0 ۳۲0۵۰ ",ماهتا مطانهه رصهمرهک رتاعامتاط مماهاهع۱11 2004 مطا عمتن فاگتامن قامصصعصه 0مامتاوه 
2009(۰) موتصامگللهن روتوم عصزه‌م منم مر معلقبامطتوظ مصتاه‌کنا وه اتمصامن 

3 ۱ متس ومامطمههه موم ۵ لامنا ربا مقلتطو۷۵ و.۲ مهتصمصاا ری ومتاعا2 و.ل مکمک .11 
1997 ), 109-121 00۰ ,37 ,۷۵۱ روطمتا ص۲0 صرح فلامو. .۰ معلقبامطانهع تعام-وتتطفییگ؟ 

(2008) وکهمصل٩‏ نصتایعظ .ممه‌عمتومه معلمیامطایمه اهءتصطعامع0 .1 10۳/212 .12 

مماممن۱1 ,16 تال ۵۶ ۲۵۵۵۲۲ ومصهوونم۳۳۵۵۵00 .0 رقتقک و٩‏ بلق و1 م1002 م5 مصهتصصه‌طمصوورهدن .13 
و12 .0 و9 ,۷۵۱, ۱9۵106011082 معلمیاوطنه۴ مه ممامصوزمگ ۵۶ امصتامل ملمتمطایهه مهومرع قامسیاوامبظ 
(2007) ,00۰73-864 

8 . م0120108 0و9 . تا ومصهونمموممم5 . فامو۱۷ . مصصم رز .. مطلهه‌مومومهن .. ,14 
۴ 080۲و ولا مه تتموری؟ رقلفنمطانهه فلت لبم 2010 عطا ۶و ومصهکعنقصه0ع۲۵۵ 
2010(۰) محططوع 1 ۵5500121100 

اجه بامط10 ۵۱۷۲۸9۰0 عصتیی میتی فصامممهصن که ممتاماجهم۲۱ ک باواقطصن0] بو مقل .15 
زممحادنا .عمزههصنوصه ملفنوطایجه وه ععصعوکصمی 0اممه 5 امه ۵۶ معهنقهعع۳۳۵ خص :2011 منم 11 ۵۶ 
(2012) 

۶ ط0۵10) ۱۵ انا مه اهتاممام۳ ممتامولمنمنا .1 مقصیم 220 ۲۰ ره/رتصطمصا لا ویک م101602" .16 
م۷ "۲۵۵1066011082 معلقیامطتو مه ومتمصهه رن [۳۵1 برقظ)ی 2 012اصامرظ: ص و ابا فصنامتم‌مهمنا 
1993(۰) ,365-379 .00 ,۳۲۵55 ]۱۸۷۲ 

۶ ۲۱۵220100 جهن فافع ماطاه1 عصللقطاو بک رفص مضه ,۲۱ رقل2۳72٩‏ .و مئژتالهاا بل 22256 .9 ,۷252 .17 
۵ ممطی/۷۷۵ هروا .۲۲.۵ طاد رفعصتلععع۱۵ظ ".عوقو لاتلمعظ ۵۶ ممتامومنامننا ما میت وعم۳ 160تباظ 
مطمنا1۱0۱0۵۵0 ٩011‏ )فصتجوه وهتاوهمم‌اصامن ممع ومتاتاه۳2 مصتاملن1 ۵ ممزدهن1 اصهاکلفی؟ معلقمنمط‌نه۴ 
,«(1995) ,]7 ,00۷/01۳0 

وه ممتمناو فصامم مگ تمزمطهها ناملا ۷۲۰ هعومک 4صه .0 ومتاعاع]۱۷ ریگ م1056 .18 
-97 .00 و1 .0 و37 ,۷۵۱ م8ط۳۵۱02009 مج ولژمو, .ملممطانعه متس اامو مصت‌صمه ۵ ممتامدگمنا110 
(1997) ,108 

0 1 ۳1۵۵ با و ۵۶ ۷۵۱۵۵ مامتا ۵ ومتاوزمام۰۵2/۵2 ۲۰ رمک 20 .۱۱ رماع ویک مققطع1 .19 
جع فطا ۶و فعصتلععع۲۳۵ ,2008 ما18 ظ مصمتاهمناطانا لهتمعمژ اهمنمهمامعن 0‏ .ص00 
معلفمتمطاه؟ لمهنصطهعا0ع0 - 2008 وومتهطمن .۰ ۲۷ . ومتجصممن۱ . زرم . لمع . مصلهعمتمص ‏ ملمنمط‌نه۴ 
.(795)2008] بضص) ,۹2012106000 ب181 ط9ت روعتصصهه 1۱ ام 4صه مصتهعصمط۳۲ 

۱۷۵0۵1۵ ممباگتتاصمن" .۱ رتمهم ومه 220 .نا رعط۷2 ویک و۱110 ویک ,100 .1 متتام۸0 ویک ,۸02116 .20 
آمع ۳ ۵۲ سول ۱۳۵۵۲۵۵ ۱.۳مصصا۲. وجابا۲ 0ععمصصحه هه ۵۶ صمنوهنا اتمتاف؟ متصصواه 10۲ 
(23-35)2003 رک660660۳0۳ 1۱ و]۷]0061( 

متصعاهو ۵۶ عصاام0 م1 لهاتم .۲ رمناا 2 بن) رقاتباظ ونژ بتالنا و.آ م2ا7/202 2 و.۲ ات۷۱0 لا معطنا .21 


سال بیست و هفتم شمارهُ دی ۱۳۹۵ نشریة مهندسی عمران فردوسی 


۱1۲ مدل‌سازی بالازدگی تونل و شریانهای... 


مه 66۵۱۵06۵ 0۲ ۵یاهل "رلزمو عاصاحصمبونا صا. مصتامموام . متمصهت م12 ۵۶ . تمزقطوه 
٩0. 12, 00. 1092-1101, )2003(‏ ,129 ۷۵۱۰ نفظ روطاهع و 6۵۵6۲۵ 

مصامرم متا ۵۶ )ممجهمهامر ورن ناوت 060ص -صمتامهگمنا من ۵۶ جمتام۳۱۵01 ویک متتامطه 1 200 .1 ر۹2821 
.(2004) 191-198 ,۲۲6۵۲5 منصصوزع مه ۷۷/۱۵۵ چم امصوظ صوورج ۲۲-۲ عصتامع۱۷۲ اصامل طا36 ۲ .عمتنتامنال 

18 عماج تبون را فعستمنای 4صامتعملصا معقا. ۵۶ مقمممروع۲ متصطفزمو رک رع۹018 تک .بلط متا 
,32 ۷۵۱۰ رکع60۱۵0) 00 کعا۵۱۵ را فصمتاهاامه مکلمنمطانعه آهمناه۷ مه آممم2تعمط مه 0ع۵ع[هاناه 
(2005) ,223-244 .00 ,4 .0 

0 160تیاها عمرژم 0۶ کتوررامصه تصمصجما ما۳ .1 مقاج‌طاه 2 20 .۲ رتتط۱۷]۵ ویاط پتالنا وبا متلاگ وگ معطلن1 
و2 ۷۵۱۰ و 160۳ ۵0 موم | آمدیمل .۱۵202 ملمنمطنهه مه امعزماناه کتعمومل 5۵ 5202060 
(2008) ,1-17 .00 و1 0۰ 

حرط عامصصی 1 «مااحطه. گم تمزمطهظ ممتگتامنا م۳1 ,5.۲ متصتمعوم۲ 4مصصمهطم)۱۷ ۱۱۲ قمع .۷۲ رن۵22 
٩۵266 ۰1600001081۷, ۷‏ صامتعع0ه۱ فصه مصتاامصصت 1 ".ععصتلوما متامین تفص واتمو. ماماهزگمنامن فط 
(2010) .158-167 ,(25)2 

ملمنمطانهه مهم اادوورمل مهو 0عاحساهه صا ومصتاه‌مرژی 0عتتباها ۵ مقصممرومم نامنا .کفتقاظ رطم۹28602720 
.(2011) ,1378-1384 .00 و10 .31,۱0 ۷۵۱۰ ,۳9۱۱66۳۱8 ع موم ۵0 یدمص اتمی ,عصن102 

و وا۵1و مااهزنن ۱1۵ جز 0مزیباها فمصتاه‌مرام ۵۶ وتورلمصه عصمصماه فنص مه مصتاممصظ مهبئتطاصمن وبا بفتاه 
2001(۰) ,۵1۷۵۲616۷ 12ماحصینامن روط ,12 

عصتاوع 1 تماه؟1۵00 ون‌ضصاط۲ ۰۷ ٩.‏ با مطات مه ,۱۷۲ ,۷۲ مطلقو۳۱۵ ویک یک رطهتهط)۱۷۲۷۲۵۱۵۵ ریک رتامصطلتتض 
فص ما کتممرم رهتصمگامت م۳ رم0۳۵۵۲۵ بوما0ط166 ۲۵۲۵۰ ۲۳6 ۵۵0۵۲۲۰۲ هاج۱ 9011 رصوتع۵:ظ 
(021100)1992عا0 ۳ ومصعلع5 اهصمتاع 

(1 9042۳ ۷۵۲5108 اب۱0۵1 00(۲8]]]]71]۷۳۴ هون ظطل۱" ررظ ,فورظ , امفطم۱۷ رجافه‌ظ 
2011(۰) ,۳۵۵۲۱۵۵۲ "9100 ۷۷۵ ۲۱۱۷۲ ۲2802 ج۵ ۱۷۲۵061 ۷۵تناتاوجمن) 

۶ صمتامانصنگ آهمتمصصا مج . ما۲ روداعمه:1. 1 ۸1ظ مطا ۵۶ مصتام۱۷۵0 ممتمممن۳ . صمصت ندال تام 
(2010) ۷۵۲616۷ ۱۵۷15 وتععط) له یام فص بگمنمان؟ ما قامصصت ]1 

لاه 660160661 ۵0۳۵۵0:۵7 .حمتاقتاآهه فلا فصع مملامع1گاقصعل وععتاو ,,۸.ظ رفظ فصه رفک ات۳۲2 
(66,)2004 181-1 :41 

مصلنهء‌متوم معمفه16 رلمیامه۱ مصتلمن و همه میمتاومن ۵۶ فلورواممه صهاع‌مه‌توم1 افد۲. له به. منت 
2001(۰)و952 م)0ومصص۱ عصتااناعممن 

معط 5 ماه مضه ۵ ممتامطصمتهنا متام مه صملامهامناوننا رهظ رق4صهتعهقعمن 
۸۵601102 روعتنظ ومصمباظ ,معا وصممم ۲۵ 8ص ععنممطمع۷ آز60 08 0۵۵۲۵۵6 ۱۱1۵۲۵0 
(1975) 

0 60و0۵ قامصصت 1 ۵۶ [ام مهبم مطا گم ولو امه ومامصصونه۳۳ ربنم رتافف؟ مه ,5.۷ رت‌صهته]۷ . 

-333 .00 و ی ۲۳2820096 ر4 .۱0 و25 ,۷۵۱ ر۵1۱9۳۵۵۱9 مساو ۳۵۵۲۵۵۵0۴8۵ ۲ رعازمگ ماماهتگهبامن] 
(2012) ,346 


22. 


23. 
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نشریهٌ مهندسی عمران فردوسی سال بیست و هفتم. شمارهدی ۱۳۹۵ 


علیرضا باقریه, ایمان لولوئی» امیرحسین باقریه ۱۳۳ 


۶ مصاام0مه لومزمصصی" و1 صنامل‌طا وبا و2ع000221 ب متفه رک رتاععظ ویک .و وکلتهط ۲.۷ م۲۳06 .35 
۰ ,41 ,۷۵۱ نامگ 0۵600606 م0۳۵0 و قاععا ممتطاصوه اه مموتنهم‌مصمی مه صمتامه۵ا110 


193-211, )2004( 


سال بیست و هفتم شمارهُ دی ۱۳۹۵ نشریة مهندسی عمران فردوسی 


۱۳۴ مدل‌سازی بالازدگی تونل و شریانهای... 


ِ فتمی شمارهٌدی ۳۹۵ 
نشریهٌ مهندسی عمران فردوسی سال بیست و هفتم شمارءٌدی ۱۳۹۵ 


